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Abstract 

We report a mechanism which can lead to long range attractions between 
like-charged spherical macroions, stemming from the existence of metastable 
ionized states. We show that the ground state of a single highly charged 
colloid plus a few excess counterions is overcharged. For the case of two 
highly charged macroions in their neutralizing divalent counterion solution we 
demonstrate that, in the regime of strong Coulomb coupling, the counterion 
clouds are very likely to be unevenly distributed, leading to one overcharged 
and one undercharged macroion. This long-living metastable configuration in 
turn leads to a long range Coulomb attraction. 
PACS numbers: 81.72Dd, 61.20Ja 
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One of the great challenges in the theory of charged colloidal suspensions is the un- 
derstanding of effective attractions between like-charged macroions that have recently been 
observed experimentally in confined systems |1|, 0], and for which no clear theoretical ex- 
planation is available. The usually employed mean field DLVO theory ^ foresees purely 
repulsive electrostatic forces between like-charged macroions. However, with divalent coun- 
terions present, simulations (using a pair of macroions) find short range attraction for high 
macroion volume fraction in aqueous systems [^] or at extremely low dielectric constant a 
Coulomb depletion force 0. Recent simulations of similar systems in aqueous solutions 
also find attractive forces [^Q. However all simulations have in common, that the ob- 
served attraction occurs only for very small distances away from the colloid surface (order 
of counterion size). 

In this letter we investigate highly charged macroions in bulk and present two important 
new results. The first concerns the ground state of an isolated macroion surrounded by excess 
counterions where it is found that the first few overcharging counterions lower considerably 
the energy. As a second finding we demonstrate that for two highly charged macroions 
separated by intermediate distances thermal fluctuations are sufficient to distribute the 
counterions unevenly, leading to one overcharged and one undercharged macroion. This 
results in a long range effective Coulomb attraction between the macroions. 

Consider one or two spherical macroions of radius and bare charge Q = —Z^e (where 
e is the elementary charge and > within the framework of the primitive model 
surrounded by an implicit solvent of relative dielectric permittivity . The small counterions 
with diameter a and charge +Zce are confined in a cubic box of length L, and the macroion(s) 
are held fixed. The colloid volume fraction is defined as A^m4vrr^/3L'^ (where is the 
number of macroions). In the case of an isolated macroion, it is located at the center of the 
box, whereas in the case of a macroion pair, they are placed symmetrically along the axis 
passing by the two centers of opposite faces. 

The (MD) method employed in the present work is similar to the one used by Kremer and 
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Grest |jTT|. To simulate a constant temperature ensemble, the ions are coupled to a heat bath 
and their motion is governed by the Langevin equation : r| = —^Vtotirt) — ""^r^ + 
fi{t), where m (chosen as unity) is the mass of the counterions, i is the i^^ counterion, 
Vfot is the total potential force made up of a Coulomb term and an excluded volume term, 
which are both pairwise additive, T is the friction coefficient, and fi a random force. These 
two last quantities are linked by the dissipation-fluctuation theorem < fi{t) ■ fj{t') >= 
QmTkBT5ij5{t — t'). For the ground state simulations the random force was set to zero. 

Excluded volume interactions are introduced via a pure short range repulsive Lennard- 
Jones (LJ) potential given by 



VLjir) = 4e 



(7 / a 



e, for r - ro < r^ut, (1) 



■.r — To/ \ r — ro , 

and otherwise, where tq = for the counterion-counterion interaction, tq = 7a for the 
macroion-counterion interaction, Vcut (= 2^^^a) is the cutoff radius. This leads to = 7.5a, 
whereas the closest center-center distance of the small ions to the macroion is therefore 
a = 8a. The Coulomb potential between a Zi and a Zj valent ion at distance r, where 
i and j denote either macroion or counterion, is given by Vcouiii^) = ksTolB^^, with the 
Bjerrum length Ib = e"^ /ineoerkBTQ, where eo is the vacuum permittivity. To link this to 
experimental units and room temperature we denote e =/cb^o ( = 298 K) and fix a = 3.57 
A. We neglect hydrodynamical interactions and hydration effects. Being interested in strong 
Coulomb coupling we choose, for the rest of this paper, = 16, corresponding to Ib = lOcr. 

To study the possibility of overcharging a single macroion, we recall the Gillespie rule also 
known as the valence-shell electron-pair repulsion (VSEPR) theory From this one knows 
that the ground state structure of two, three and four electrons disposed on a hard sphere 
corresponds to simple geometrical situations, namely a line (two electrons diametrically 
opposed), a triangle, and a tetrahedron, respectively. A straightforward calculation shows 
that, for a central charge of +2e , the maximally obtainable overcharging is -2e (i.e. 2 
electrons), being independent of macroion radius. The excess electrons gain more energy 
by assuming a topological favorable configuration than by escaping to infinity, the simple 



reason of overcharge. We resort to simulations to elucidate this behavior for one colloid with 
a high central charge. 

To quantify this phenomenon, we have considered three macroionic charge Zm'- 50, 90 
and 180 corresponding to a surface charge density of one elementary charge per 180, 100 and 
50 A^, respectively, and fixed Z^. = 1 for the rest of this letter. Then we add successively 
overcharging counterions (OC). The electrostatic energy as a function of the number of OC 
is displayed Fig. 1. We note that the maximal (critical) acceptance of OC (4, 6 and 8) 
increases with the macroionic charge (50, 90 and 180 respectively). Furthermore for a given 
number of OC, the gain in energy is always increasing with Zm- Also, for a given macroionic 
charge, the gain in energy between two successive overcharged states is decreasing with 
the number of OC. Note that at T = 0, the value acts only as a prefactor. It means 
that the ground state structure is solely dictated by topological rules (i.e., the counterions 
arrangement around the sphere). 

The resulting curve can be very simply explained by assuming that the energy e per ion 
on the surface of a neutralized macroion depends linearly on the inverse distance between 
them, hence is proportional to ^/N for fixed area, where N is the number of counterions 
on the surface. The energy gain AEi = {N + l)e{N + 1) — Ne{N) of the first OC is a 
pure surface correlation term. For the next OC one needs to take into account the Coulomb 
repulsion IsZ^/a, leading to lowest order in 1/A^ for the energy gain of the n*^ OC: 



AEn = ne{N) 
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Determining e(A^) from the measured value for Ai^i, we obtain a curve that matches the 
simulation data almost perfectly, compare Fig.l. 

An energy per ion, which scales like \fN , has been found for an ionic Wigner crystal 
(WC) on a planar surface where each ion interacts with an oppositely charged background 
charge which is smeared out over its Wigner-Seitz cell. This energy per ion is given by 
e{c)/kBTo = —ac^^^lBZ"^, with a = 1.96, and c is the two-dimensional concentration of the 



crystallized counterions of valence Z^. [|T^]. This Ansatz has been tried recently to explain 



strong ionic correlations observed in various soft matter system • In our simulation we 



find for ^Ei/ {kBTo) -18.0, -24.4, and -35.3 for = 50, 90, and 180, respectively, whereas 
the Wigner crystal scenario predicts -21.0, -28.0, and -39.5, which is off by a decreasing rate 
of 17 — 12%. This might be due to the assumption of a homogeneous background charge, 
and the assumption of a planar geometry, neither of which are fully fulfilled, however the 
error gets smaller for higher values of Zm- 

Using the Wigner crystal ionic energy and eq. (§), the maximally obtainable number 
^mocc of OC counterions is readily found to be 
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This value depends only on the number of counterions A^. It originates from the topological 
arrangement of the ions around a central charge, and is independent of Bjerrum length or 



radius of the macro ion. For large Q it reduces to the form Q*maxl^ ~ * \/ZmZc which 



was derived in Ref. |]T5[ in a more elaborate fashion. 

To obtain the interaction potential profile, we added one counterion coming from infinity 
towards a macroion of bare charge Zm = 180 and computed the global electrostatic energy of 
the system, see Fig. 2. The first OC starts to gain correlational energy at distance r 12a 
from the center of the colloid, which is about Aa from the surface. This fits only roughly 
with the distance Z^Ib/A predicted from WC theory |jl^,|l5|, and is more of the order c~^/^. 



With adding more excess counterions the Coulomb barrier increases, and for the ninth OC 
it exceeds the gain in correlational energy, when being on the macroion surface. Thus the 
configuration becomes metastable. The curve for the first OC can be nicely fitted with an 
exponential fit of the form Ei{r) /kBT^ = — 35.3exp [— 7.1(r — a)/a]. For the n^^ OC simply 
the appropriate Coulomb monopole contribution AlB{n — l)/r needs to be added, see Fig. 2. 
This exponential dependence is not predicted by the WC theory, where a 1/r dependence 
should be seen due to the interaction of the removed ion with its correlation hole. 

Next, we consider two spherical like charged macro ions at a colloidal volume fraction 
fm = 7-10~^ at room temperature Tq, at fixed center-center separation R, in presence of their 



divalent counterions (ensuring global charge neutrality). Initially the counterions are ran- 
domly generated. Figure 3 shows two macroions surrounded by their quasi-two-dimensional 
counterions layer. The striking peculiarity in this configuration is that it corresponds to an 
overcharged and an undercharged sphere. There is one counterion more on the left sphere 
and one less on the right sphere compared to the bare colloid charge. Such a configuration is 
referred as ionized state. In a total of 10 typical runs, we observe this phenomenon 5 times. 
We have also carefully checked against a situation with periodic boundary conditions, yield- 
ing identical results. However it is clear that such a state is in "pseudo-equilibrium" because 
it is not the lowest energy state. 

To estimate the energy barrier, electrostatic energy profiles at zero temperature were 
computed, where we move one counterion from the overcharged macroion to the under- 
charged, restoring the neutral state (see Fig. 4). We have checked that the path leading 
to the lowest barrier of such a process corresponds to the line joining the two macroions 
centers. One clearly observes a barrier, which increases linearly with the charge The 
ground state corresponds as expected to the neutral state. The overcharged state is only 
slightly higher in energy, the difference being approximately the monopole contribution 
E/ksTQ — Ib{4:/S — 4/12) 1.67. The physical origin of this barrier can be understood 
from the single macroion case where we showed that a counterion gains high correlational 
energy near the surface. This gain is roughly equal for both macroion surfaces and decreases 
rapidly with increasing distance from the surfaces, leading to the energy barrier with its 
maximum near the midpoint. For the single macroion case we showed that the correlational 
energy gain scales with V^nj whereas here we observe a linear behavior of the barrier height 
with Zm- We attribute this effect to additional ionic correlations since both macroions are 
close enough for their surface ions to interact strongly. For large separations we find again 
that the barrier height increases with V%^, as expected. This dependence of the barrier 
also shows that at room temperature such ionized states only can occur for large Z^- In our 
case only for = 180, the ionized state was stable for all accessible computation times. 
Unfortunately, it is not possible to get a satisfactory accuracy of the energy jumps at non- 
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zero temperatures. Nevertheless, since we are interested in the strong Coulomb coupling 
regime, which is energy dominated, the zero temperature analysis is sufficient to capture the 
essential physics. 

Results concerning the effective forces at zero temperature between the two macroions 
are now investigated which expression is given by 

Pe.ff{R) = Fmm{R) + FlJ + Fmc , (4) 

where Fmm.{R) is the direct Coulomb force between macroions, Flj is the excluded volume 
force between a given macroion and its surrounding counterions and Fmc is the Coulomb 
force between a given macroion and all the counterions. Because of symmetry, we focus 
on one macroion. To understand the extra-attraction effect of these ionized-like states, we 
consider three cases: (i) Fjo„ = Ff.ff in the ionized state (ii) F^eut = -^e// in the neutral case 
(iii) Fmono = Feff simply from the effective monopole contribution. Our results are displayed 
in Fig. 5 for Z.^ = 180, where the ionized state was also observed at room temperature. 
The non-compensated case leads to a very important extra attraction. This becomes drastic 
for the charge asymmetry of ± 2 counterions at short separation R/a ^ 2.5, a situation 
which was also observed in our simulation at room temperature [|rB|. In contrast to previous 



studies [|, ^, these attractions are long range. For a sufficiently large macroion separation 
(from 3.5) the effective force approaches in good approximation the monopole contribution. 

In summary, we have shown that a sufficiently charged colloid can in principle be highly 
overcharged due to correlation effects of the counterions, and that this effect is quantitatively 
well described by a Wigner crystal, i.e. eqs. (0) and (H). In the strong Coulomb coupling 
regime, this energy gain can be of the order of many /cbTq. 

Furthermore, due to this energetically favorable overcharged state it was found that for 
two like-charged macroions, an initially randomly placed counterion cloud of their neutral- 
izing divalent counterions may not be equally distributed after relaxation, leading to two 
macroions of opposite net charges. The resulting configuration is metastable, however sep- 
arated by an energy barrier of several /cbTq when the bare charge is sufficiently large. Such 



configuration possess a natural strong long range attraction. 

We acknowledge useful discussions with H. Lwen and M. Deserno. This work is supported 
by Laboratoires Europeens Associes (LEA) and a computer time grant hkf06 from NIC 
Jiilich. 
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List of Figures 

FIG. 1: Electrostatic energy (in units of /cbTq) for zero temperature configurations of a 
single charged macroion of radius = 7.5(7 as a function of the number of overcharging 
counterions for three different bare charges Q (in units of e ) . The neutral case was chosen 
as the potential energy origin, and the curves were produced using the theory of eq.(2), 
compare text. 

FIG. 2: Electrostatic energy (in units of /cbTq) of a divalent counterion as function of 
distance from the center of a macroion with radius = l.ba and charge Q = —180 (in 
units of e ). The energy is normalized to zero at distance infinity. Data and fits are shown 
for the first, the second, the eighth and ninth overcharging (OC) counterion. 

FIG. 3 Snapshot of a "pseudo-equilibrium" configuration at room temperature Tq where 
the counterion-layers do not exactly compensate the macroions charge. Here the deficiency 
charge is ±1 counterion (or ± 2e as indicated above the macroions) and R/ a = 3.6. 

FIG. 4 Total electrostatic energy (in units of /csTq) of the system, for zero temperature 
configurations, of two macroions at a center-center separation ofi?/a = 2.4asa function of 
one displaced counterion distance from the left macroion for three typical values Q (in units 
of e ). The exact neutral state was chosen as the potential energy origin. The lines are a 
guide to the eye. The insert indicates the path (dotted fine) of the moved counterion. The 
ending arrows of the arc indicate the start position (left sphere) and final position (right 
sphere) of the moved counterion. 

FIG. 5 Reduced effective force between the two spherical macroions at zero temperature 
for Zm — 180 as a function of distance from the center. The different forces are explained 
in the text. The lines are a guide to the eye. 
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